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ABSTRACT 

An algebraic turbulence model is proposed for use 
with three-dimensional Navier-Stokes analyses. It in- 
corporates features of both the Baldwin- Lomax and 
Cebeci-Smith models. The Baldwin- Lomax model uses 
the maximum of a function f(y) to determine length 
and velocity scales. An analysis of the Baldwin-Lomax 
model shows that f(y ) can have a spurious maximum 
close to the wall, causing numerical problems and non- 
physical results. The proposed model uses integral re- 
lations to determine 6*u e and 6 used in the Cebeci- 
Smith model. It eliminates a constant in the Baldwin- 
Lomax model and determines the two remaining con- 
stants by comparison to the Cebeci-Smith formula- 
tion. Pressure gradient effects, a new wake model, 
and the implementation of these features in a three- 
dimensional Navier-Stokes code are also described. Re- 
sult are shown for a flat plate boundary layer, an an- 
nular turbine cascade, and endwall heat transfer in 
a linear turbine cascade. The heat transfer results 
agree well with experimental data which shows large 
variations in endwall Stanton number contours with 
Reynolds number. 

INTRODUCTION 

The Baldwin-Lomax turbulence model (BLM) was 
developed in 1981 [1] as a simple algebraic turbulence 
model for use with Navier-Stokes codes. The model is 
based on the widely-used Cebeci-Smith boundary layer 
model (CSM, [2]), with modifications made to avoid 
the need to search for the edge of the boundary layer. 

Since its development, the BLM has become one 
of the most widely used turbulence models in com- 
putational fluid dynamics. NASA’s RECON library 
database lists 287 references written since the original 
publication [1] that included the names Baldwin or Lo- 
max in their abstracts. Most references indicated that 
the model was successful. Some noted physical details 
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missing from the model, e.g. pressure gradient effects 
[3] which are discussed later. A few noted disagree- 
ment between the constants in the BLM and the CSM 
[4, 5, 6]. Finally, several references noted numerical 
problems with the model, usually a failure to find the 
correct outer length scale. This problem is covered in 
detail here and an alternative model is proposed. 

The present work began as an attempt to calcu- 
late endwall heat transfer in a linear turbine cascade 
studied experimentally by Boyle and Russell [7], who 
found two distinctly different heat transfer patterns at 
low and high Reynolds numbers. Calculations of the 
cascade were made with the three-dimensional Navier- 
Stokes code described in [8, 9]. Initial calculations us- 
ing the BLM looked promising but showed isolated re- 
gions of unrealistically high or low heat transfer. The 
problem was traced to multiple peaks in the Baldwin- 
Lomax function f(y) used to determine a turbulent 
length scale, which led to incorrect estimates of the 
turbulent viscosity and boundary layer thickness. 

An analysis of the BLM was performed, similar to 
[4, 5, 6]. Like the previous work, it was found that the 
constants used in the BLM do not give agreement with 
the CSM. Unlike the previous work, it was also shown 
that f(y) can have a spurious maxima at the edge of 
the viscous sublayer, especially in favorable pressure 
gradients and at low Reynolds numbers. 

A new turbulence model was developed that avoids 
searching for / mar . It uses the relation / 0 fdy = 6*u e 
(shown independently in [10]) directly in the Cebeci- 
Smith model and correlates $ with the centroid of the 
f(y) curve. The model avoids one of the constants 
used in the BLM, and determines the two remaining 
constants by comparison with the CSM. It also includes 
pressure gradient effects and a new wake model. 

Calculations made of the endwall heat transfer prob- 
lem using the new model show the change in heat trans- 
fer pattern with Reynolds number and agree fairly well 
with the experiment. 

NUMERICAL METHOD 

The code described by Chima in [8, 9] was used in 
the present work, and is described below briefly. The 
code, referred to as rvc3d (Rotor Viscous Code 3-D) 
is a Navier-Stokes code for analysis of steady flows in 



turbomachinery blade passages. The governing equa- 
tions are formulated in Cartesian coordinates including 
rotation about the x-axis, and transformed to general 
body-fitted coordinates. The Navier-Stokes equations 
are discretized using second-order central finite differ- 
ences. The thin-layer approximation is used to neglect 
streamwise viscous terms, but all blade- to- blade and 
hub-to-tip viscous terms are retained. 

An explicit four-stage Runge-Kutta time-marching 
scheme is used to solve the flow equations. Fourth- 
difference artificial dissipation terms are added. A spa- 
tially varying time step and implicit residual smoothing 
are used for convergence acceleration. 

C-type grids are generated in a blade-to-blade plane 
using a version of Steger and Sorenson's GRAPE code. 
The C-grids give good resolution of the leading edge 
and the wake. Grids are stacked spanwise using a hy- 
perbolic tangent stretching function. 

Boundary conditions for the cascade heat transfer 
problem were as follows: On the blade and endwall 
the no-slip condition was imposed on the velocity com- 
ponents, the wall pressure was found from the normal 
momentum equation, and the wall temperature was set 
to 1.2 x Toin to cause heat transfer from the walls. The 
flow was assumed to be periodic blade-to-blade and 
symmetric about the mid-span. At the exit boundary 
p was specified and /?, pu , pv , and pw were extrapolated 
from the interior. 

' At the inlet boundary, Pq , To, v and w profiles were 
specified. The velocity profiles were computed using a 
turbulent law- of- the- wall profile matched to the mea- 
sured inlet boundary layer thickness and freestream ve- 
locity. P 0 was computed from the velocity profile and 
the measured static pressure. To was computed us- 
ing a thermal law-of-the-wall profile chosen to match 
the measured inlet Stanton number. Stanton numbers 
within the blade row were insensitive to the inlet value. 
A Riemann invariant R~ based on u is extrapolated 
from an interior grid point to the inlet. A correction 
to R~~ is made for non-isentropic flow due to bound- 
ary layer dissipation and endwall heat transfer, using 
a technique to be described in a later paper. R~ and 
To are then solved for u and c at the inlet. 

CEBECI- SMITH MODEL 

In Navier-Stokes analyses the effective viscosity and 
effective thermal conductivity are taken to be the sum 
of laminar and turbulent components, i.e., 

Pe// — P/am T Pturb (1) 

P/am = / T 

Pre/ \T re j ) 

L n* ( P/<*m Pturb ^ 

*•" = c ' ( -p r +_ sr ) 

where Pr f = 0.9 is the turbulent Prandtl number and 
n = 0.667. 


The Cebeci-Smith model [2] is a two-layer algebraic 
model for the turbulent viscosity: 


Pturb — 


{ 


Hi for y <y c 
Ho for y > y c 


( 2 ) 


where y c is the first point at which Hi > Po- 

A Prandtl-van Driest formulation is used for the in- 
ner layer: 
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Here p is the density, k = 0.4 is the von Karman 
constant, D is the van Driest damping function with 
A + = 26., and \du/dy\ is the absolute value of the 
velocity gradient. 

A Clauser formulation is used for the outer layer: 


Po — Kpj6*u e 


(4) 


where 6* is the displacement thickness and u e is the 
edge velocity. The Klebanoff intermittency function is 
given by: 
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BALDWIN-LOMAX MODEL 

The Baldwin-Lomax Model [1] is derived from the 
Cebeci-Smith model described above. In the BLM the 
inner formulation is identical to (3) except that the 
magnitude of the vorticity |u>| replaces \du/dy\. 

The BLM replaces the outer formulation of the CSM 
with: 


Po = KpjC cp min 


{ 


Vmax Jmax 

^wkVmax j j / fm ax 


(6) 


where K = .0168 is the Clauser constant and C cp — 1.6 
is an additional constant. 

The product y ma xfmax replaces 6*xt e in the CSM. 
The function f{y) is defined by: 


m 

= y\u\D 

(7) 

fmax 

= max(/(y)) 


Vmax 

— * 2/ 3-f fmax 



Later it is shown that y ma x usually occurs at a fixed 
fraction of 6, but that in certain cases y m ax occurs at 
the edge of the viscous sublayer, and the BLM fails. 

In the intermittency function (5) the BLM replaces 
y/S with Cxieb y/y m ax> where Cxieb = 0.3. In wakes 
D is set to 1 and the lower option in (6) is used (it 
may also be active at separation.) C w k was originally 
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reported as 0.25 but is generally taken to be 1.0. The 
term Udijj = |Kn ar | — | is the difference between 
the maximum and minimum total velocities in a profile. 
Recent papers have replaced V ma x with V(y ma x )■ 


BALDWIN-LOMAX MODEL ANALYSIS 

Several authors have analyzed the BLM to show why 
the function f(y) is used, and to adjust Caleb and C cp 
for agreement with the CSM [4, 5, 6]. The analysis 
of [5] is given briefly below for reference. It is then 
extended to show when the BLM might fail, and fi- 
nally used to develop a new model based on the original 
CSM. 

Cole’s turbulent velocity profile is used in the anal- 
ysis of the BLM. In the viscous sublayer: 

u + = y + (8) 


and in the wall- wake region: 
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U e 


— In [ Re^ 
K 


n 
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( 9 ) 


where B = 5.5, and II is a pressure gradient parameter 
given by: 


II =0.9(/?+ .5) 75 , /? = _l_dp 1 (10 ) 

7~u>all dX 

An infinitely favorable pressure gradient has II = 0, 
an equilibrium gradient has II = 0.5, and an infinitely 
adverse gradient has II = oo. For a specified value of 
II, (9) may be evaluated at u = u e and y = <5, and 
solved for the friction velocity u*. The equation has 
two roots with opposite signs, giving separated profiles 


when u* is negative. Figure 1 shows several turbulent 
velocity profiles over a range of II for Res = 105,000, 
the value used in [5], 

The Baldwin-Lomax function f(y) = y|u/|D can be 
computed from (9) by replacing |u>| with \du/dy\ and 
neglecting D, which is meant to damp spurious maxima 
in the viscous sublayer (y/S < 0.05). The resulting 
expression for f(y) is: 

/W = »H = M[, + n^.i„(f)] on 

The maximum of f(y) occurs where | sin (^) is max- 
imum, which occurs independently of II at 

ymax — 0.6466 (12) 

Figure 2 shows profiles of the Baldwin-Lomax func- 
tion f(y) for the velocity profiles shown in figure 1. 

For an infinitely adverse gradient (II = oo) T wa ii — 
y+ = D = fi Q — 0, which can cause numerical prob- 
lems. For an adverse gradient (e.g. II =5) there is a 
clear maximum in f(y) and the BLM is well-behaved. 
For an equilibrium layer (II = 0.5) there is a weak max- 
imum in the outer part of the boundary layer where 
grids are relatively coarse and small errors in the solu- 
tion can lead to large errors in t/ ma x- For an infinitely 
favorable gradient (II — 0) f(y) is constant in the outer 
layer and the BLM fails. 

The constants Caleb and C cp are evaluated by com- 
paring the BLM to the CSM. Briefly, Caleb is eval- 
uated by comparing the intermittency formulations, 
giving Caleb = 0.646. Baldwin and Lomax suggested 
Caleb = 0.3, giving little intermittency. C cp is eval- 
uated by comparing the outer formulations. The dis- 
placement thickness 6* is found by integrating (u e — u) 
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using (9). The result is: 


Ccp .646(1 + 1.819711) ^ 

Thus C cp is a strong function of pressure gradient, 
ranging from C cp = 1.548 at II = 0 to C cp = 0.851 at 
II = oo. The equilibrium value C cp = 1.216 at II = 0.5 
is recommended over Baldwin and Lomax’s suggested 
value of C cp = 1.6. 


In practice r is replaced by fi \u\. Results using the two 
formulations are discussed later. 

Cebeci’s extension of van Driest ’s damping model [2] 
is used to allow the boundary layer to relaminarize in 
response to favorable pressure gradients. 

i 26 , v dp 

~ s/l + 11.8p+ ’ p ~ p^ds 

The local pressure gradient is evaluated using: 


SPURIOUS MAXIMA IN f(y) 

For the BLM to work properly, y ma x must reliably 
correlate with <5, as in equation (12). Several authors 
have reported problems with spurious maxima in f(y) 
very close to the wall, e.g., [4, 10]. Here it is shown that 
under some conditions f(y) has a spurious maximum 
in the viscous sublayer. 

The function f(y) given in (11) is based on the wall- 
wake profile(9). The function can also be evaluated in 
the viscous sublayer using (8), giving: 

D = y + D (14) 


ftu 


sub 


du + 


dy+ 


dp 

ds 



(18) 


The “edge” velocity V e is evaluated at a grid line far 
enough from the wall to indicate the general flow di- 
rection. 

The outer formulation is given by: 


Po 
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Kpj min 


1+5.5 


/ F 

( C w ky{\^max 
^ CKleb 6 


\-\Vmin\) 


(19) 


the maximum of which occurs at the edge of the sub- 
layer. 

For simplicity consider an infinitely favorable pres- 
sure gradient (II = 0). The sublayer and wall- wake pro- 
files meet at = 11.635, where (14) gives f+ ax — 
4.2. The outer profile (11) gives /+ ar = 1/k — 2.5. 
Thus for II = 0, f(y) has two local maxima, and 
the largest occurs at the edge of the viscous sublayer, 
even when the van Driest damping term D is used to 
suppress the spurious maximum. In reality a buffer 
layer blends the viscous sublayer and the wall region 
smoothly, and the spurious maximum is reduced. Sim- 
ilar results are found for certain II > 0, depending on 
Re*. For example, for a fairly high Re* = 105,000 
(fig. 1 and 2) f max — sub ^ f max — outer for II <C 0.28. It 
can be shown that spurious maxima are most likely in 
favorable pressure gradients and at low Re*. 

PROPOSED TURBULENCE MODEL 

The inner formulation is similar to the BLM formu- 
lation: 

m = pi 2 M (15) 

where /, D y and are as defined in (3) and k = 0.4. 

The friction velocity u* used in (3) is usually defined 
in terms of the wall shear r^au. Several researchers 
have modified u* by replacing T wail = 0 with the local 
shear in the profile, either to avoid problems at separa- 
tion [4] where r wa u = 0, or as an indirect way to allow 
the boundary layer to respond to pressure gradients 

[ii]- 



where CKieb — 0.55 and C w k = 0.825, as derived below. 

The function F is an integral of the BL function f(y) 
without the damping term D. 

f°° du 

F = Jo J V ~dy dy ~ ^ Ue ( 20 ) 

shown using integration by parts and assuming that 
limy_**(<9u/dy) = 0. This relationship was discovered 
independently by Johnson [10] but was rejected in favor 
of another scheme based on f max . In practice du/dy 
is replaced with |u;|, which goes to zero at S , even if 
du/dy does not. The integral scale F provides an un- 
ambiguous measure of the product wit hout having 
to evaluate either term directly. It is used in the outer 
formulation of the CSM (4), avoiding the constant C rp 
entirely. 

The length scale y is the centroid of the f(y) curve, 
defined by: 



f(y)dy 


( 21 ) 


Substituting (11) for /(y), the integrals may be evalu 
ated analytically. The result is a transcendental equa 
tion for y/S as a weak function of n, which is solved 
numerically. For II = 0, y/S = 0.5, which is ob- 
vious from figure 2. For n = .5, y/S — .55, and 
for n = oo, y/S — .606. The equilibrium value 
y/S — .55 = Cxieb was used in the present work. 

The model is implemented by replacing du/dy with 
|w| in (20) and integrating numerically. The accumu- 
lated integral is stored in an array at each point in the 
profile. When the integration is complete, Ff 2 gives 
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FLAT PLATE VELOCITY PROFILES 
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Figure 3. — Comparison of flat plate velocity profiles to 
Spalding’s composite law of the wall. 
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Figure 4. — Computational grid for the annular turbine 
cascade. 


the value F at y, and a search through the array gives 

m- 

The wake constant C wk is evaluated by equating the 
two options in (19) assuming a separating profile (IT = 
oo), giving: 


C wk = 


Fsep 

VsepAV 


( 22 ) 


Substituting F aep = u e 6/ 2, y sep = .6066, and A V/u e « 
1 (valid for separation or in the near wake) gives C wk = 
0.825. 

In this work AV was taken to be — iVTnml, as 

implied in the original BLM. It may be better to use 
AV = IV’(y)!- \V m i n \ , as in recent work with the BLM. 

Finally, the inner and outer turbulent viscosities are 
combined using Granville’s blending function [12]: 


He/f = Ho tanh 


Hi 

Ho 


(23) 


3-D IMPLEMENTATION 

The turbulence model is applied in the blade-to- 
blade (tj) and spanwise (£) directions independently. 
The length scale y is computed differently in the inner 
and outer layers. For the inner layer the Buleev length 
scale based on the distance from the hub and blade is 
used to model near-wall effects. 


Vi 


2s 


n s C 


+ s < + 


+ s$ 


(24) 


Using the local shear modification, the inner layer be- 
comes independent of profile direction. 

For the outer layer the actual distance from the ap- 
propriate wall is used to model boundary layer thick- 
ness effects, i.e., y 0 = s n or s The outer layer gives 


different profiles of Hturb in the rj and £ directions. The 
two profiles are combined vectorally: 

Hturb — yHtr) + Ht{ ( 25 ) 

Although there is no physical argument for this com- 
bination, it has the appropriate behavior at both the 
wall and the free stream. 

RESULTS 

Calculations were made of a two-dimensional flat 
plate boundary layer at Moo = 0.25. The grid had 248 
points on the plate and 95 points normal to the plate, 
with the second grid point at y+ < 1. Figure 3 shows 
three velocity profiles at Re x = 1.5 x 10 6 compared 
to Spalding’s composite expression for the law of the 
wall. Profiles computed with both the BLM and the 
present model using u* = y/r wa u/p agree well with 
Spalding’s expression. However the profile computed 
with the present model using the local shear modifi- 
cation u* = ^/r(y)/p does not agree with the law of 
the wall. Apparently this disagreement has not been 
noticed before. It is recommended that the local shear 
modification be used with caution. 

An annular cascade of turbine vanes tested experi- 
mentally at NASA Lewis Research Center by Goldman 
and McLallin [13] was analyzed. The computational 
grid had 127 x 445 x 63 = 360, 045 points and is shown 
in figure 4. The solution converged 3.5 decades in 1,800 
iterations, which took about 2.4 hours on a Cray Y-MP. 

The vane has a constant cross section but the dis- 
tance between vanes increases with radius, thereby in- 
creasing the vane loading. This is seen in the surface 
pressure distributions shown in figure 5. The computed 
pressures agree very well with the experimental data. 
The spanwise variation of mixed out total pressure loss 
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STATIC PRESSURE DISTRIBUTION 



FRACTION OF AXIAL CHORD 

Figure 5. — Computed and measured pressure distribu- 
tions for the annular turbine cascade. 



Figure 6. — Computed and measured loss coefficient 
profiles for the annular turbine cascade. 


coefficient downstream of the vanes is shown in figure 6. 
The calculated losses are slightly higher than the mea- 
sured losses. Figure 7 compares calculated and mea- 
sured efficiency contours across the wake at a distance 
of 1/3 axial chord downstream of the trailing edge. The 
kinetic energy efficiency is defined by: 

V = 2C p (7o - T) (26) 

where Q is the total velocity, To is the inlet total tem- 
perature, and T is the static temperature. The over- 
all measured efficiency was 0.96 and the computed ef- 
ficiency was 0.94. The shape of the computed wake 
is very similar to the measured wake, showing under- 
turned, high loss regions near the endwalls due to sec- 
ondary flows. However the computed wake has not 
mixed out as much as the measured wake, suggesting 
that the wake turbulence model may need improve- 
ment. 

Calculations have also been made of endwall heat 
transfer in a linear turbine cascade. The cascade was 
tested experimentally by Boyle and Russell [7] over a 
wide range of chord Reynolds numbers Re c and inlet 
boundary layer thicknesses <5; n . Endwall heat trans- 
fer was determined by appling electrical power to a 
foil heater covered by a liquid crystal sheet. The tem- 
perature at which the liquid crystal sheet exhibited a 
specific color was known from a calibration test. Col- 
ored bands on the sheet indicated isotherms, and since 
the heat flux was uniform they were also lines of con- 
stant heat transfer. By photographing the isotherms 
at several power settings the endwall Stanton number 
pattern was mapped. The Stanton number is defined 
by: 

- (k^L) 

Cf __ V dn/wall (*27} 

~ ( pV)i„C p (T wall - T oin ) K ] 


It was found that St increased somewhat as <5,- n in- 
creased , with little change in the overall contour pat- 
tern. However, two distinct St contour patterns were 
seen at low and high Re c , as discussed below. 

Two heat transfer cases were computed using the 
new turbulence model. The cases had similar inlet 
boundary layer thicknesses but very different Re c . The 
local shear modification to the van Driest damping (16) 
was used for both cases before the disagreement with 
the law of the wall was noted. The effects of this mod- 
ification are discussed at the end of this section. The 
computational grid used for both cases is shown in fig- 
ure 8. It has 165 x 44 x 41 = 297,660 points on the 
half span. The second grid point was at < 0.7 for 
the low Re c case and at < 3.5 for the high Re r case. 

The first case had Re c = 78,000, S in = .2 x span, 
and Mi n ~ 0.03. To improve the convergence rate 
of the code, Mi n was scaled to 0.095 while maintain- 
ing Re c . The code took 4,000 iterations to converge 
three decades. Total CPU time on a Cray Y-MP was 
4.4 hours. Figure 9 compares measured and computed 
endwall St contours for the low Re c case. The exper- 
imental St contours (right) increase through the pas- 
sage as the flow accelerates. The contours run from the 
leading edge region of the upper blade along the pres- 
sure surface, then deflect gradually towards the suction 
surface, resembling secondary flovc streamlines near the 
endwall. The computed St contours (left) are about 20 
percent low at the inlet but show fairly good agreement 
with the measured patterns in the front of the passage. 
The computed contours do not follow the pressure sur- 
face as closely as the measured contours, but the peak 
St in the throat region agrees reasonably well with the 
data. 

The second case had Re c = 490,000, <5 ln = .15 x 
span, and = 0.175. The code converged about 3.5 
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Figure 7. ^Computed and measured contours of kinetic energy efficiency in the wake of the annular turbine 
cascade, turbine cascade. 


orders of magnitude in 2.8 CPU hours. Figure 10 com- 
pares measured and computed endwall St contours for 
the high Re c case. Experimental St contours (right) 
are lower than the values seen at low Re c and follow a 
different pattern. Here the experimental contours run 
almost directly from blade to blade, and the computed 
contours (left) agree closely. Endwall heat transfer aug- 
mentation near the blade leading edge is well predicted. 

The two cases were analyzed to determine the 
physics behind the different St contour patterns. It 
was found that the differences are probably not due to 
differences in either the presure field or in secondary 
flow patterns, which are very similar for the two cases. 
The endwall Stanton numbers were quite sensitive to 
modifications to the van Driest damping. With a con- 
stant value of A+ = 26 the Stanton numbers were high 
but did show some differences between the St patterns 
at different Re c . The local shear modification (16) 
reduced Stanton numbers by about 25 percent every- 
where for both cases, giving better quantitative agree- 
ment with the experimental data. Both the local shear 
modification (16) and Cebeci’s modification (18) accen- 
tuated the slope of the contours towards the front of 
the passage in the low Re c case but had little effect at 
high Re c . The differences in the Stanton number pat- 
terns are thus consistent with relaminarization of the 
near-wall flow at low Reynolds numbers. It is hoped 
that further modifications to the van Driest damping 
formulation will better model the effects of favorable 
pressure gradients while remaining consistent with the 
law of the wail. 


SUMMARY 

A three-dimensional Navier-Stokes code was used to 
predict endwall heat transfer in a turbine cascade. Ini- 
tial calculations using the Baldwin-Lomax turbulence 
model showed erratic heat transfer patterns due to 
a spurious maximum in the Baldwin-Lomax function 
f(y) in the viscous sublayer. An analysis of the model 
showed that the problem was most likely to occur at 
the low Reynolds numbers and favorable pressure gra- 
dients found in low speed turbine cascades. 

A new turbulence model was developed that uses 
unambiguous integral relations to estimate 6*u e and 
<5. The model includes modifications to the van Driest 
damping function to account for the effects of pressure 
gradient. A new wake model was also proposed. Model 
constants were determined by comparison with the 
Cebeci-Smith model. It should be simple to retrofit the 
new model into existing codes that use the Baldwin- 
Lomax model. 

For a flat plate boundary layer calculation, both the 
Baldwin-Lomax model and the present model showed 
good agreement with Spalding’s expression for the law 
of the wall. However a commonly-used modification 
to van Driest’s damping function caused disagreement 
with the law of the wall. 

Calculations of flow in an annular turbine cascade 
showed good agreement with experimental pressure 
distributions. Computed wakes looked similar to mea- 
sured wakes but showed less mixing, suggesting that 
the wake model may need improvement. 

Calculations were also made of endwall heat transfer 
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in a linear turbine cascade at two Reynolds numbers. 
Comparisons with experimental data showed that the 
model predicts observed variations in end wall Stan- 
ton number patterns with Reynolds number reasonably 
well. 
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